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1.  INTRODUCTION 


\\ 

J  Generalized  networks  are  an  Important  class  of  optimization 
models,  with  uses  in  a  wide  variety  of  fields.  This  report  describes 
the  development  and  implementation  of  a  generalized  network 
algorithm. 

In  133,  Jarvis  et.  al.  recommend  a  generalized  network  model, 
system  for  closure  optimization  and  planning  (SCOPE),  for  crisis 
action  deployment  planning.  In  SCOPE,  large  generalized  networks 
must  be  repeatedly  solved.  These  networks  have  special  structure, 
which  results  in  computational  advantages. 

In  this  report,  a  generalized  network  implementation  is  developed 
for  solving  very  large  generalized  networks.  This  Implementation 
Includes  new  data  structures  for  storing  the  basis,  in-core/out-of - 
core  handling  of  the  arcs,  and  special  handling  of  pure  network 
structure. 

In  this  report,  a  detailed  examination  of  the  SCOPE  model  is 
provided  and  its  ef f ect  on  implementation  Issues  is  d ' scussed.  The 
SCOPE  model  is  highly  structured.  This  report  demonstrates  how  this 
structure  can  be  used  to  advantage.  In  a  companion  report .[43, 

- - - .  .  ^  j? 


extensive  testing  is  presented  which  addresses  the  question  *What 
affects  the  computation  time  for  a  SCOPE  model 7”^ 

Section  2  provides  an  Introduction  to  modeling  with  pure  and 
generalized  networks.  Section  3  gives  an  overview  of  linear 
programming,  as  it  applies  to  generalized  network  solution  methods. 

In  section  4,  the  special  structure  of  a  generalized  network  basis  is 
detailed  and  efficient  storage  methods  are  developed.  Section  3 
gives  methods  for  handling  the  arc  data.  This  includes  both  storage 
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Methods  and  tachnlquas  for  determining  the  arc  to  enter  the  basis.  A 
sethod  for  determining  the  arc  to  leave  the  basis  is  given  in  section 
6.  Section  7  describes,  in  detail,  the  algorithms  used  to  update  the 
structures.  Section  8  discusses  some  related  implementation 
concerns.  These  include  the  initial  basis  to  be  used,  and  the  effect 
of  embedded  pure  network  structure.  Section  9  gives  a  review  of  the 
conclusions. 
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2.0  MODELING  WITH  GENERALIZED  NETWORKS 


Generalized  networks,  aa  the  name  suggests,  are  a  generalization 
of  standard,  or  "pure",  networks.  By  modifying  a  restriction  that 
occurs  in  pure  networks,  many  previously  Intractable  problems  have 
been  modeled.  Section  2. 1  provides  a  brief  overview  of  pure  network 
modeling.  Section  2. 2  expands  this  overview  to  include  generalized 
networks.  Section  2. 3  details  the  SCOPE  generalized  network  model. 
The  SCOPE  model  is  an  example  of  generalized  network  model.  It  will 
form  the  basis  for  the  examples  used  later  in  this  report. 

2. 1  Pure  Networks 

A  "pure*  network  can  be  thought  of  as  a  pipeline  system.  This 
system  has  suppliers  and  users  of  the  materiel  flowing  through  the 
pipe.  Each  of  the  suppliers  and  users  has  a  known  supply  and  demand 
respectively.  The  pipeline  connects  the  suppliers  and  users.  The 
pipeline  may  have  intermediate  junction  points  which  are  not 
suppliers  or  users.  (An  example  would  be  a  pumping  station. > 

Each  pipe  has  a  known  capacity.  This  capacity  represents  the 
limit  on  the  pipe  expressed  in  the  rate  of  flow  of  materiel  through 
the  pipe.  There  is  a  unit  cost  associated  with  materiel  that  flows 
through  the  pipe.  This  cost  is  linear  in  the  quantity  of  materiel) 
in  other  words,  if  the  amount  flowing  through  a  pipe  is  doubled,  then 
the  cost  is  also  doubled. 

The  objective  is  to  move  the  materiel  from  the  suppliers  to  the 
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ui*ra  through  the  pipeline  at  minimum  cost.  This  involves  assigning 
flows  to  pipes  so  that  (1)  each  user  gets  the  aaount  needed;  (2)  no 
supplier  sends  more  than  available;  (3)  no  pipe  has  sore  materiel 
than  its  capacity;  (4)  no  materiel  enters  or  leaves  the  system  except 
at  users  or  suppliers;  and  (3)  vhat  enters  the  pipe  at  one  end, 
leaves  it  at  the  other. 

The  suppliers,  users,  and  junctions  where  two  or  more  pipes  come 
together  may  be  represented  by  points,  called  nodes.  The  lengths  of 
pipes  between  points  (nodes)  are  called  arcs.  Associated  with  each 
node  is  a  number,  called  its  requirement.  If  the  requirement  is 
negative,  then  the  node  is  a  supplier,  and  the  number  is  the  amount 
that  it  can  supply.  If  the  requirement  is  positive,  then  the  node  is 
a  user,  and  the  number  represents  the  amount  it  demands.  A  zero 
requirement  can  be  used  for  junctions  where  arcs  meet  without 
representing  a  supplier  or  user. 

Associated  with  each  node  is  a  flag.  The  flag  indicates  whether 
the  requirement  must  be  met  exactly,  or  whether  the  absolute  value  of 
the  requirement  represents  an  upper  bound  for  the  supplier  (or  user). 
This  accommodates  models  with  nodes  in  which  suppliers  must  ship  the 
full  amount  of  their  supply  and  for  users  that  may  or  may  not  use  the 
full  demand  Indicated  for  them. 

Arcs  have  capacities,  cost  and  direction.  The  node  that  an  arc 
begins  at  is  called  its  "TAIL”.  The  ending  node  is  the  "HEAD".  For 
arc  number  ARC  these  two  ends  are  referred  to  as  TAIL (ARC)  and 
HEAD (ARC)  respectively. 

Figure  2-1  shows  a  sample  network.  The  circles  represent  nodes. 
The  lines  (arrows)  between  them  are  the  arcs. 
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2.2  Generalized  Hetvorka 


Generalized  networks  are  extensions  of  pure  networks.  The 
difference  is  that  the  restriction  that  "flow  into  an  arc  equals  flow 
out  of  the  arc"  is  relaxed.  Instead,  a  generalized  network  allows 
for  "leaky”  arcs  or  arcs  that  gain  flow.  The  loss  or  gain  is 
specific  to  each  arc  and  can  vary  throughout  the  network.  The  only 
requirement  is  that  the  arc  must  gain  or  lose  a  constant  fraction  of 
its  flow.  For  instance,  an  arc  might  always  triple  its  flow;  another 
pipe  might  always  quarter  its  flow.  The  fraction  that  the  arc 
changes  the  flow  is  called  the  multiplier.  The  multipliers  for  the 
example  arcs  are  3  and  0. 25  respectively.  As  a  matter  of  convention, 
costs  are  calculated  by  multiplying  the  flow  that  enters  the  arc  by 
the  arc  cost. 

Because  of  gains  and  losses  in  a  generalized  network,  it  is 
impossible  to  require  that  supply  equals  demand,  as  is  normal  in  pure 
network  models.  "Slack"  and  "surplus”  arcs  are  required  to  model  the 
excess  supply  or  demand  that  normally  occurs  in  a  generalized 
networks.  These  "slack"  and  "surplus"  arcs  are  modeled  using  "self- 
loops"  at  each  node.  "Self-loops"  are  arcs  that  begin  and  end  at  the 
same  node.  Surplus  arcs,  associated  with  demand  nodes,  have 
multipliers  of  *1.  Slack  arcs,  attached  to  supply  nodes  have 
multipliers  of  -1. 

All  other  properties  of  pure  networks  hold  for  generalized 
networks.  A  generalized  network  with  all  multipliers  equal  to  one  is 
a  pure  network.  An  arc  with  a  multiplier  of  one  is  called  a  pure 
arc. 

Generalized  network  models  are  very  useful  in  many  problems. 
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1)  In  a  financial  model,  materiel  flowing  in  a  network  repreaenta 
aoney  and  nodea  repreaent  varioua  pointa  in  time.  Hultipliera  can  be 
used  to  represent  increase  in  money  due  to  interest. 

2)  In  an  energy  allocation  model,  materiel  flowing  in  a  network 
represents  electricity,  and  arcs  represent  physical  wires. 

Multipliers  can  be  used  to  model  energy  losses  that  result  when 
electricity  flows  along  a  wire. 

3)  In  a  deployment  model,  materiel  flowing  in  the  network  are 
men  and  equipment  to  be  moved.  Arcs  represent  movement,  by  either 
airplane  or  ship.  If  materiel  is  moved  by  air  then  the  weight 
(STONS)  of  the  materiel  determines  the  amount  that  can  be  moved.  If 
movement  is  by  sea,  then  the  volume  (HTONS)  or  square  footage  of  the 
movement  requirement  is  critical.  Multipliers  can  be  used  to  model 
conversion  of  weight  into  volume. 

2.3  The  SCOPE  MRMATE  Model 

In  C33,  a  method  for  solving  a  large  deployment  problem  was 
presented.  This  method,  called  System  for  Closure  Planning  and 
Evaluation  < SCOPE),  addressed  the  following  problem: 

Given  a  set  of  assets  (airplanes  and  ships),  a  set  of  movement 
requirements  (people,  ammunition,  etc.),  and  a  set  of  ports  to  use, 
is  there  a  way  to  move  the  movement  requirements  with  the  available 
assets  through  ports  so  that  the  requirements  arrive  at  the  target 
area  when  needed? 

A  much  more  detailed  examination  of  the  problem  is  provided  in 
[33.  The  method  proposed  is  based  on  three  optimization  components: 
a  network  flow  with  side  constraints  which  assigns  assets  to  pairs  of 
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porta;  •  generalized  network  flow  which  aaalgna  the  novenent 
requirements  to  assets;  and  a  Benders'  constraint  generator  to  link 
together  the  two  models.  This  report  is  concerned  only  with  the 
second  of  the  optimization  pieces. 

Fundamental  to  the  SCOPE  method  is  the  concept  of  channels.  A 
channel  consists  of  a  pair  of  ports,  together  with  a  number  of 
Identical  assets.  One  port,  the  Port  of  Embarkation  (POE),  is  where 
the  movement  requirements  will  be  loaded  onto  the  assets.  The  other 
port,  the  Port  of  Debarkation  (POD),  is  the  destination  of  the 
assets.  The  assets  are  assumed  to  cycle  between  the  POE  and  POD. 
Depending  on  the  distance  between  the  POE  and  POD,  the  assets  may  be 
able  to  make  one  or  more  trips  between  the  POE  and  POD  in  a  single 
time  period  (which  can  be  taken  as  a  day  for  simplicity).  Or,  if  the 
ports  are  far  apart  or  the  asset  moves  slowly,  it  may  take  several 
days  to  cycle  between  the  POE  and  the  POD.  The  capability  of  the 
channel  is  the  rate  at  which  the  assets  deliver  movement  requirements 
to  the  POD  from  the  POE. 

The  first  optimization  model  in  SCOPE  determines  the  channel 
capabilities.  The  second,  generalized  network  model,  must  assign  the 
movement  requirements  to  channels  at  specific  time  periods  so  as  to 
meet  strategic  objectives. 

Suppliers  in  the  SCOPE  model  are  the  movement  requirements. 

Users  are  time  expanded  channels.  Time  expansion  refers  to  the  fact 
that  a  channel  can  move  a  given  amount  on  each  day.  Nodes  will  be 


created  for  the  channel  for  each  day. 

Certain  movement  requirements  cannot  be  moved  on  specific 
channels.  For  instance,  a  movement  requirement  may  not  be  air 


transportable,  so  they  cannot  use  channels  using  air  assets.  There 
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are  many  other  restrictions.  There  is  ■  certain  delay  Involved  in 
getting  a  Movement  requirement  to  the  POE,  so  even  if  a  movement 
requirement  can  use  a  channel,  it  may  not  be  able  to  use  it  on 
certain  days.  Arcs  are  created  from  each  movement  requirement  to 
each  time  expanded  channel  node  that  the  movement  requirement  can 
use.  (See  Figure  2-2). 

The  interpretation  of  flow  on  an  arc  from  movement  requirement  M 
to  channel  C  on  day  D  is  that  amount  of  M  vill  arrive  at  the  POE 
associated  vlth  C  on  day  D  and  vill  use  the  assets  associated  with  C 
to  be  transported  to  the  POD  of  C.  From  this  information,  it  is 
possible  to  determine  vhen  M  vill  arrive  at  its  final  destination. 

The  cost  associated  vith  an  arc  depends  on  the  value  of  getting  M  to 
its  final  destination  at  that  time. 

There  are  many  mays  of  assigning  costa  to  the  arcs,  but  one  of 
the  most  flexible  involves  time  vlndovs.  Each  movement  requirement 
has  a  vindov  of  days  in  vhich  it  is  desired  to  arrive  at  its  final 
destination.  If  it  arrives  at  its  final  destination  vithin  its 
vindov  there  is  no  cost.  If  it  arrives  outside  its  vindov  (either 
early  or  late)  then  the  cost  is  a  function  of  hov  many  days  early  or 
late  it  arrives. 

The  final  complication  is  in  hov  to  measure  the  size  of  a  flov. 
There  are  tvo  types  of  channels:  air  and  sea.  The  major  limitation 
on  the  amount  an  air  channel  can  move  is  veight  (STOHS)  of  the 
movement  requirements.  The  major  limitation  on  the  sea  channels  is 
volume  (HTQNS).  Each  movement  requirement  has  a  veight  and  volume, 
and  the  relationship  betveen  these  tvo  values  depends  on  the  movement 
requirement. 

To  model  this,  the  veight  of  the  movement  requirement  is  taken  as 


1  can  use 

2  can  use 

3  can  use 

4  can  use 


Channels  1 
Channels  2 
Channel  3 


and  2 
and  3 


•SCOPE  Basic  Model 


Its  size.  If  the  movement  requirement  is  eent  by  an  air  channel,  no 
converalon  takes  place.  If  the  movement  requirement  la  sent  by  a  sea 
channel,  the  weight  is  converted  to  volume  by  a  multiplier  on  the 

arc.  For  instance,  if  a  movement  requirement  weighs  20  STONS  and  has 

volume  SO  flTONS,  arcs  to  air  channels  will  be  pure  arcs  (have 

multiplier  one)  and  arcs  to  sea  channels  will  have  multipliers  of  2.5 

(>50/20).  Figure  2-3  provides  the  complete  network  for  the  example. 

The  network  created,  called  the  MRMATE  network,  has  significant 
structure.  The  major  features  aret 

1)  all  nodes  are  either  suppliers  or  users.  In  fact,  this  network 
is  a  transportation  network  (see  Bazaraa  and  Jarvis  C 1 J ) . 

2)  many  of  the  multipliers  are  one. 

3)  all  arcs  out  of  a  movement  requirement  have  one  of  two 
multipliers. 

3)  many  arcs  have  zero  cost. 

4)  the  arcs  have  no  capacities. 

The  largest  problem  that  could  occur  in  practice  is  estimated  to 
have  approximately  2000  nodes  and  over  500, 000  arcs.  Furthermore, 
these  problems  must  be  solved  repeatedly  in  a  very  short  amount  of 
time. 

The  MRMATE  model  will  be  used  as  an  example  of  the  type  of 
specializations  possible  in  implementations  of  generalized  networks 
in  the  remainder  of  this  report. 
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3.0  SOLVING  GENERALIZED  NETWORKS 

Sine*  the  generalized  network  model  la  ao  uaeful,  it  la  Important 
to  have  a  computer  aolutlon  technique  that  will  find  aolutiona 
quickly.  It  can  be  ahovn  that  generalized  netvorke  are  a  apeclal 
caae  of  linear  programming,  ao  that  any  technique  to  aolve  linear 
programming,  like  the  Simplex  method,  can  be  uaed  to  find  aolutiona 
to  thia  model.  But  there  are  dlaadvantagea  of  uaing  theae  general 
purpoae  algorithma.  The  beat  known  methoda  take  too  much  apace  and 
are  relatively  alow.  Fortunately,  the  Simplex  method  can  be 
apeciallzed  ao  aa  to  take  advantage  of  the  apeclal  atructure  In  a 
generalized  network.  The  apeciallzed  aimplex  method  aolvea 
generalized  networka  quickly  uaing  little  apace. 

Section  3. 1  Linear  Programming 

Linear  programming  la  the  moat  fundamental  optimization  model  in 
operationa  reaearch.  Bazaraa  and  Jarvia  (13  provide  an  excellent 
Introduction  to  thia  field.  The  linear  programming  model  employe  the 
optimization  of  a  linear  function  eubject  to  a  aet  of  linear 
conatrainta.  A  linear  function  la  a  function  that  la  of  the  form: 

C»  x,  ♦  C,  x.  ♦ .  ♦  C„  x. 

where  each  of  C( ,  C« ,  ...  C*  are  conetanta 
and  x( ,  x« ,  and  x,  are  the  varlablee. 


A  linear  conatralnt  la  of  the  form: 

A,  X,  ♦  Aa  xa  ♦ .  ♦  A,  X.  ■  B 

where  A, ,  Aa ,  ...  A.  and  B  are  all  conatanta. 

When  a  generalized  network  ia  repreaented  aa  a  linear  program. 
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variable  x,  ia  associated  with  each  arc.  This  variable  represents 
the  amount  of  t lov  in  the  arc. 

There  is  a  linear  constraint  for  each  node.  This  constraint 
controls  the  amount  of  the  flow  that  exits  or  enters  the  node.  The  B 
value  for  the  constraint  is  the  supply  or  demand  for  that  node.  The 
constraint  forces  the  net  amount  of  flow  at  a  node  (including  the 
self-loop),  to  be  the  supply  or  demand  for  the  node. 

Every  arc  la  associated  with  just  two  nodes:  the  tail  and  head 
nodes  for  the  arc.  This  implies  that  the  variable  associated  vith 
each  arc  occurs  in  at  moat  two  constraints  in  the  linear  program. 
Self-loops  occur  in  only  one  constraint. 

In  matrix  terms,  each  constraint  represents  a  row  of  the 
constraint  matrix  and  each  variable  represents  a  column.  The 
preceding  argument  indicates  that  there  are  at  most  tvo  non-zero 
elements  in  each  column  of  the  constraint  matrix. 

Because  the  definition  of  generalized  networks  in  this  report 
allows  Just  a  single  multiplier  for  each  arc,  one  of  the  non-zeros  of 
each  column  will  be  the  multiplier  on  the  arc.  For  arcs  that  are  not 
self-loops,  the  other  non-zero  element  of  the  column  will  be  -1.  The 
multiplier  will  be  in  the  row  associated  with  the  constraint  on  the 
head  of  the  arc.  The  -1,  for  non-self -loops,  is  associated  with  the 
tall  of  the  arc. 

Every  variable  is  assumed  to  be  constrained  to  be  nonnegative. 
There  are  standard  "tricks"  to  transform  variables  not  of  this  form 
to  the  assumed  form. 

Since  the  arcs  (variables)  have  capacities,  it  is  necessary  to 
treat  this  problem  as  a  linear  program  with  upper  bounds.  These 
upper  bounds  are  linear  constraints  themselves.  Due  to  the  simplicity 
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of  upper  bound*  It  la  possible  to  treat  thes  implicitly  in  the 
solution  algorithm. 

3. 2  Solving  Linear  Programs  -  The  Primal  Simplex  Method 

There  are  many  methods  for  solving  linear  programs.  The  most 
widely  used  is  the  primal  simplex  method.  This  technique  has  proved 
to  be  efficient,  both  in  execution  time  and  computer  space. 

For  every  linear  program,  there  is  an  optimal  solution  with  no 
more  than  one  non-zero  variable  for  each  constraint.  This  is 
referred  to  as  a  basic  optimal  solution.  The  optimal  non-zero 
variables  form  a  basis.  A  basis  is  any  set  of  variables  with  the 
following  properties: 

1)  There  are  not  more  variables  than  constraints  in  the  linear 
program. 

2)  No  column  of  the  constraint  matrix  for  any  variable  in  the 
basis  can  be  written  as  a  weighted  cbmbinatlon  of  the  columns  of  the 
other  variables  in  the  basis. 

3)  There  is  a  feasible  solution  to  the  linear  program  using  Just 
the  variables  in  the  basis. 

The  steps  of  the  primal  simplex  method  are  as  follows: 

1)  Find  an  initial  basis. 

2)  Find  a  variable,  not  in  the  basis,  to  enter  the  basis.  If  none 
exists,  STOP.  The  current  basis  is  optimal. 

3)  Find  the  variable  in  the  basis  that  will  leave. 

4)  Update  the  basis 

3)  GOTO  Step  2. 

One  iteration  of  steps  2  through  4  is  called  a  "pivot”. 
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Step  1  can  ba  accompliahed  In  varioua  vaya.  Tha  aiaplaat  aathod 
la  to  taka  tha  *alack"  and  *aurplua*  varlablaa  that  oftan  occur  in  a 
linaar  program  and  uaa  them  aa  tha  Initial  baaia.  Somatimaa 
artificial  varlablaa  auat  ba  addad  vhara  "real*  alack  and  aurplua 
varlablaa  do  not  axiat.  Thaaa  artificial  varlablaa  ara  givan  a  high 
coat,  ao  that  tha  optimal  aolutlon  will  not  amploy  any  of  tham. 

It  la  oftan  poaaibla  to  datarmina  a  aat  of  varlablaa  that  craataa 
a  vary  good  aolutlon.  Thia  uaually  raducaa  tha  numbar  of  pivota 
raquirad  to  raach  optimality.  Tha  tlma  to  find  a  good  atartlng 
aolutlon,  called  an  advanced  atart,  auat  ba  ahort  enough  not  to 
offaat  tha  reduced  computation  time  for  tha  raat  of  tha  algorithm. 

Identifying  a  variable  to  enter  tha  baaia  ia  accompliahad  by 
determining  tha  change  in  objective  function  if  tha  variable  ia 
lncraaaad  by  a  amall  value.  Thia  change  la  called  tha  reduced  coat. 

Incraaaing  tha  value  of  tha  variable  entering  tha  baaia  vill 
change  tha  valuaa  of  tha  currant  baaic  varlablaa.  One  of  thaaa 
varlablaa  vill  be  tha  firat  to  raach  zero.  Thia  la  tha  variable  to 
exit  tha  baaia. 

Tha  nov  baaia  conalata  of  tha  old  baaia,  without  tha  exiting 
variable,  and  tha  entering  variable.  Varioua  valuaa  muat  ba  updated, 
including  tha  nav  variable  valuaa  and  the  reduced  coata  for  varlablaa 
not  in  tha  baaia. 

Tha  primal  aimplax  method  can  ba  adapted  for  upper  bounda  on 
varlablaa.  Rather  than  treat  tha  upper  bounda  aa  "normal* 
conatralnta,  which  would  ba  Inefficient,  tha  definition  of  baaia  la 
alightly  redefined.  A  non-baaic  variable  can  nov  have  value  of 
either  zero  or  ita  upper  bound.  A  baaic  variable  can  have  any  value 
between  zero  and  ita  upper  bound.  A  non-baaic  variable  at  ita  upper 
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bound  way  tnttr  tha  basis  if  dscrsasing  its  flow  slightly  improves 
^bs  objective  function*  When  the  basic  variables  change  value  (Step 
i  3*  on®  Of  them  will  reach  its  upper  bound  or  zero  first*  That 


variable  will  be  the  exiting  variable. 


4.0  STORING  THE  BASIS 


The  main  reason  that  a  specialized  simplex  method  is  faster  than 
a  general  purpose  simplex  method  for  generalized  networks  is  that  the 
basis  has  a  special  structure.  This  structure  makes  every  simplex 
computation  easier.  This  section  defines  the  basis  structure  and 
gives  data  structures  to  efficiently  store  it. 

4. 1  Basis  Definition 

A  basis  in  linear  programming  consists  of  a  set  of  columns,  one 
for  each  rov,  with  the  property  that  no  column  1s  a  weighted  sum  of 
the  others.  In  a  generalized  network,  columns  correspond  to  arcs,  so 
the  basis  is  a  set  of  arcs.  Rows  correspond  to  nodes,  so  there  is 
one  arc  in  the  basis  for  each  node.  The  final  property,  called 
linear  Independence,  is  more  complicated  to  describe. 

If  a  set  of  arcs  is  examined,  the  set  of  nodes  will  be 
partitioned  into  seta  of  nodes  that  are  connected  to  each  other  (see 
Figure  «. 1).  These  sets  are  called  components.  Within  a  component, 
the  arcs  can  form  cycles.  Self-loops  are  treated  as  cycles  of  length 
1.  A  component  can  have  zero,  one,  or  more  than  one  cycle 
(components  A,  B,  and  C  respectively  in  Figure  4.1).  It  is  possible 
to  show  that  if  a  component  has  more  than  one  cycle,  there  is  at 
least  one  arc  that  is  the  weighted  sum  of  the  other  arcs  in  the 
component.  Therefore,  for  a  set  of  arcs  to  form  a  basis,  it  is 
necessary  that  no  component  formed  by  the  arcs  have  more  than  one 
cycle.  It  is  also  possible  to  show  that  if  a  component  has  no  cycle. 


some  other  component  must  have  more  than  one  cycle.  Therefore,  it  is 
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•Iso  ntceaoary  that  every  component  hava  at  least  ona  eye la.  So 
avary  component  has  axactly  ona  cycla.  A  coaponant  with  axactly  one 
cycle  is  called  a  "one-tree*. 

One  further  condition  la  required  to  ensure  linear  independence. 
If  the  cycle  is  not  a  self-loop,  it  is  possible  for  an  arc  in  the 
cycle  to  be  a  weighted  sum  of  the  other  arcs  in  the  cycle.  A 
necessary  and  sufficient  condition  for  this  not  to  occur  is  for  the 
cycla  to  have  a  cycle  multiplier  not  equal  to  one.  The  cycla 
multiplier  is  calculated  as  follows:  Assign  an  orientation  to  the 
cycle  (clockwise  or  counter-clockwise).  The  cycle  multiplier  is  the 
product  of  the  arc  multipliers  for  those  arcs  pointed  in  the  same 
direction  as  the  orientation,  divided  by  the  product  of  the  arc 
multipliers  of  those  arcs  pointed  in  the  reverse  direction  as  the 
orientation. 

To  summarize,  a  sat  of  arcs  is  a  basis  if  the  following 
conditions  are  satisfied: 

1 )  The  number  of  arcs  is  equal  to  the  number  of  nodes. 

2)  Each  component  has  exactly  one  cycle. 

3)  Each  component  with  a  cycle  that  is  not  a  self -loop  has  a 
cycle  multiplier  that  is  not  equal  to  one. 

Some  valid  and  invalid  basis  examples  are  given  in  Figure  4-2. 

Since  every  multiplier  in  a  pure  network  is  1,  it  is  not  possible 
for  a  pure  network  basis  to  have  a  cycle  that  is  not  a  self  loop. 

4. 1. 1  Specialization  for  the  SCOPE  Model 

The  arc  multipliers  in  the  MRMATE  model  have  a  special  structure. 
Since  every  arc  connects  a  movement  requirement  to  either  an  air  or  a 
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(a)  valid 


label:  multiplier 


(b)  valid 


(c)  invalid  (no  cycle) 


(d)  invalid  (two  cycles) 


figure  4-2 .  Valid  and  Invalid  Basis 
21 


label  on  arc:  multiplier 

(e)  invalid:  cycle  multiplier  of  1 


figure  U-2,  Continued 
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sea  channel,  the  area  have  eultipllera  of  either  1  or  the  converaion 
factor  for  the  Movement  requirement.  Given  the  vay  that  cyclea  form 
In  the  MRtlATE  model.  It  la  eaay  to  ahov  that  If  a  component  haa  just 
air  channel  area  or  just  aea  channel  area,  then  the  cycle  aeeoclated 
with  the  component  muat  be  a  aelf  loop.  In  other  vorda.  If  a 
component  haa  only  one  type  of  arc  then  vhat  component  haa  the  aame 
baala  atructure  aa  a  pure  network  baala.  Since  pure  netvorkB  can  be 
aolved  more  efficiently  than  generalized  netvorka.  It  la  likely  that 
aome  advantage  can  be  taken  of  the  baala  atructure  In  thla  caae. 

4. 2  Storing  the  Baala 

Since  the  baala  for  generalized  netvorka  haa  apeclal  atructure. 

It  ahould  be  poaaible  to  atore  the  baala  in  an  efficient  vay.  There 
are  two  Important  factora  In  atoring  the  baala s  atorage  apace  and 
computation  time. 

The  following  aectlona  outline  a  method  of  atoring  the  baala, 
called  the  linked  rooted  tree  method.  Thla  method  la  aimilar  to  that 
of  Brown  and  McBride  C23,  but  diff era  in  aome  important  vaya. 

The  linked  rooted  tree  method  la  baaed  on  the  data  atructurea 
uaed  for  pure  netvorka  (see  Xennington  ES3>.  In  thla  method,  nodes  on 
the  cycle  are  seen  as  roots  for  trees  consisting  of  nodes  not  on  the 
cycle.  These  trees  are  then  linked  around  the  cycle  (see  Figure  4- 
3).  Each  component  contains  one  or  more  trees  together  with  the 
linking  cycle. 


4. 2. 1  Predecessor  Structure 


The  moat  fundamental  operation  required  for  the  manipulation  la 
to  "go  up"  the  tree.  (Here  the  cycle  ia  conaidered  "on  top"  of  the 
tree).  Thia  ia  required  in  determining  the  arc  to  exit  the  baaia, 
for  the  area  that  muat  be  checked  are  exactly  thoae  above  the 
endpointa  of  the  entering  arc.  It  ia  alao  very  uaeful  vhen 
determining  the  nev  baaia  (Section  7). 

For  any  aelected  node  not  on  the  cycle,  the  predeceaaor  (PRED)  ia 
defined  to  be  the  (unique)  node,  auch  that  there  ia  a  baaic  arc 
connecting  the  tvo  nodea  and  PRED  ia  cloaer  to  the  cycle  than  the 
aelected  node. 

For  nodea  on  the  cycle,  an  arbitrary  orientation  of  the  cycle  ia 
aelected.  The  PRED  of  a  node  on  the  cycle  la  the  node  Juat  before  it 
on  the  cycle  uaing  the  aelected  orientation.  The  PRED  of  a  node  on 
the  cycle  la  alao  on  the  cycle.  If  the  cycle  la  a  aelf-loop,  the 
PRED  of  the  cycle  node  ia  defined  to  be  ltaelf.  (See  Figure  4-4). 

The  PRED  provldea  the  only  method  of  moving  from  one  cycle  node 
to  another  under  the  linked  rooted  tree  ayatem.  Alao  note  that  there 
la  no  connectlona  between  componenta,  for  that  ability  ia  not 
required  for  the  aimplex  calculatlona. 

4. 2. 2  Thread  Structure 
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The  thread  atructure  (THREAD)  provldea  a  mechaniam  for  vlaiting 
every  node  in  a  tree.  The  order  in  which  the  nodea  are  vlaited  la 
defined  to  be  the  "preorder  traveraal"  (aee  Kennington  [53).  Thia 
order  haa  the  property  that  if  node  X  la  on  the  path  from  Y  to  the 
root  then  node  X  ia  vlaited  before  node  Y.  Note  that  the  THREAD  ia 
only  within  treea,  not  between  them.  (See  Figure  4-5). 
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The  THREAD  la  required  by  the  beeie  update  routinee  to  deteraine 
thoae  nodes  whose  duals  and  LEVELS  (Section  4.2.3)  aust  be  updated. 

It  is  in  this  structure  that  the  linked  rooted  tree  aethod 
differs  froa  that  used  by  Brown  and  McBride.  In  that  re port,  the 
THREAD  was  defined  traverse  around  the  cycle  in  the  opposite 
direction  of  PRED. 

4. 2. 3  Level  Structure 

The  level  structure  (LEVEL)  gives  the  "distance11  of  a  node  to  the 
cycle.  Nodes  on  the  cycle  have  a  LEVEL  value  of  zero.  This 
structure  is  required  by  the  routine  to  find  the  exiting  arc  (Section 
6).  (See  Figure  4-6). 

4. 2. 4  Reverse  Thread  Structure 

The  reverse  thread  structure  (RTHREAD)  is  siaply  the  inverse  of 
the  THREAD  structure.  This  permits  the  visiting  of  nodes  in  reverse 
order.  Typically,  this  structure  is  only  used  to  aake  the  basis 
update  aore  efficient.  (See  Figure  4-7) 

4. 2. 3  Arc  Inforaation  Structures 

Some  information  on  the  basic  arcs  is  required  to  perform  the 
simplex  calculations.  This  Includes  the  CAPACITY,  arc  multiplier 
(MULT)  and  current  FLOW,  to  determine  the  arc  to  exit  the  basis;  the 
arc  COST,  to  update  the  dual  variables;  and  the  arc  number  (ARC)  to 
record  the  optimal  solution.  If  the  Inforaation  on  all  arcs  (basic 
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Reverse  Thread 


Figure  4-7.  Reverse  Thread  Structure 


and  non-baaic)  la  avallabla,  than  it  la  only  nacaaaary  to  atora  tha 
ARC  valua  axplicltly.  Aa  Sactlon  3. 1  vill  ahov,  hovavar,  for  larga 
problaaa  It  la  nacaaaary  to  hava  only  a  liaitad  aaount  of  arc 
inforaation  availabla  at  any  givan  tiaa.  Tharafora,  all  of  tha  abova 
baalc  arc  inforaation  auat  ba  atorad. 

Tha  inforaation  on  tha  baaic  arc  that  connacta  NODE  and 
PRED(NODE)  la  aaaociatad  with  NODE.  Sinca  it  la  not  claar  vhathar 
tha  arc  bagina  at  NODE  and  anda  at  PRED(NODE)  or  tha  ravaraa,  tha  ARC 
valua  ia  givan  a  aign  dapandlng  on  tha  oriantation.  In  tha  foraar 
caaa,  ARC  la  poaltlvai  in  tha  lattar,  ARC  ia  nagatlva.  (Saa  Pigura  4- 
6). 

4. 2. 6  Dual  Valua  Structura 

Tha  dual  variabla  (DUAL)  for  aach  noda  ia  raquirad  to  dataraina 
an  arc  to  antar  tha  baaia.  Tharafora,  tha  dual  for  aach  noda  ia 
ratalnad  at  all  tlaaa. 

4. 2. 7  Cycla  Hultipliar 

Tha  cycla  aultlpllar  (CHULT)  la  daflnad  aa  follovas  glva  an 
oriantation  to  tha  cycla j  tha  cycla  aultlpllar  ia  tha  product  of  tha 
arc  aultlpllara  for  area  in  tha  aaaa  dlractlon  aa  tha  oriantation 
dlvldad  by  tha  arc  aultlpllara  of  thoaa  area  in  tha  ravaraa 
dlractlon. 

CHULT  ia  daflnad  only  for  nodaa  on  cyclaa  that  ara  not  aalf- 
loopa.  CHULT  ia  tha  aaaa  for  all  nodaa  on  tha  aaaa  cycla. 

Tha  cycla  aultlpllar  can  ba  thought  of  aa  tha  aaount  of  flov  that 


label  on  arc:  arc  number 
direction  on  arc:  predecessor 


Figure  4-8.  Arc  Structure 
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will  result  if  on*  unit  of  flov  ia  sent  around  tha  cycle  in  the  given 
orientation.  In  a  valid  baeia,  CMULT  cannot  be  on*. 

CHULT  ia  needed  to  deteraine  the  arc  to  exit  the  baeia,  and  to 
update  the  flova.  Becauae  of  .the  number  of  multiplication*  and 
diviaiona  required,  thie  la  a  very  time-coneuming  number  to 
calculate.  Fortunately,  thia  number  muat  be  calculated  juat  once  for 
each  new  cycle  created.  So  pivota  that  do  not  create  a  nev  cycle  do 
not  require  the  calculation  of  any  CHULT  value*. 


5.  Handling  Area 


In  general  there  are  far  aore  area  than  nodea  in  a  network  aodel. 
The  handling  of  the  area  la  therefore  critical  to  deereaaing 
coaputatlon  tiae  and  atorage  requireaenta.  Deapite  the  large  aaount 
of  are  data,  very  little  data  la  required  for  any  individual  pivot. 
Thla  peraita  the  atorage  of  the  arc  data  outaide  central  coaputer 
aeaory,  normally  on  a  high  apeed  aaaa  atorage  device. 

Hethoda  for  handling  the  are  data  vaa  reported  in  detail  in  C  84- 
093.  The  following  aectiona  review  the  concluaiona  of  that  report. 

S. 1  Chooaing  an  Entering  Arc 

The  priaal  simplex  method  providea  flexibility  in  the  choice  of 
arc  to  enter  the  baaia.  The  only  property  the  entering  arc  auat  have 
la  that  placing  a  aaall  aaount  of  flow  on  the  arc  will  decreaae  the 
objective  function  value.  It  la  not  neceaeary  to  aelect  the  arc  that 
will  yield  the  greateat  decreaae  in  the  objective  value  the  aoat. 

Given  the  current  dual  valuea,  it  ia  eaay  to  deteraine  the  effect 
of  Baking  a  aaall  change  in  the  current  flow  of  a  non-baaic  arc.  If 
the  arc  currently  haa  no  flow  on  it,  increasing  the  flow  by  one  unit 
would  add  the  following  aaount  to  the  objective  function! 

HULT  •  DUAL (HEAD)  -  DUAL(TAIL)  -  COST 
where  the  arc  in  queation  la  froa  node  TAIL  to  node  HEAD  and  haa 
aultiplier  HULT  and  coat  COST. 

If  the  current  flow  on  the  arc  ia  the  capacity  of  the  arc,  then 
deereaaing  the  flow  by  one  unit  would  add  the  following  aaount  to  the 
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objectives 

•  (MULT  •  DUAL (HEAD)  -  DUAL (TAIL)  -  COST) 

It  eight  not  be  possible  to  change  the  flows  by  one  unit.  Some 
arc  in  the  basis  could  reach  one  of  its  bounds  before  one  unit  of 
flov  is  placed  on  the  entering  arc.  Conversely,  it  sight  be  possible 
to  change  the  flov  by  more  than  one  unit.  The  actual  amount  of 
change  is  determined  in  Step  3  of  the  primal  simplex  algorithm. 

The  above  equations  are  referred  to  as  the  reduced  cost  for  the 
arc.  If  the  reduced  cost  is  negative,  the  objective  value  would  be 
decreased  if  the  non-baslc  arc  were  to  enter  the  basis.  Non-basic 
arcs  with  negative  reduced  cost  are  eligible  to  enter  the  basis.  If 
there  is  no  arc  with  a  negative  reduced  cost  then  the  current 
solution  is  optimal. 

Generally,  there  are  many  arcs  that  are  eligible  to  enter  the 

basis  at  each  iteration.  It  is  necessary  to  choose  from  among  those 

possible.  For  instance,  it  is  possible  to  calculate  all  of  the 

reduced  costa  and  choose  the  arc  with  the  most  negative  reduced  cost 

% 

to  enter  the  basis.  This,  generally,  will  have  fever  pivots  than 
other  methods  but  the  amount  of  time  required  to  calculate  all  of  the 
reduced  costs  would  be  prohibitive. 

Another  alternative  is  to  calculate  the  reduced  costs,  one  at  a 
time,  and  choose  an  arc  to  enter  the  basis  as  soon  as  one  is  found 
with  a  negative  reduced  cost.  This  will  minimize  the  amount  of  time 
to  calculate  reduced  costs i  but  it  will  cause  many  pivots  which 
improve,  only  marginally,  the  objective  value. 

In  practice,  two  methods  are  used  to  choose  an  arc  to  enter  the 
basis.  These  methods  are  referred  to  as  the  fixed  page  method  and  the 
candidate  list  method. 
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5.1.1  Fixed  Pag*  Method 


In  the  Fixed  Page  method,  the  reduced  coat  lor  a  fixed  number  of 
arcs  (a  page)  la  calculated  at  each  pivot.  The  arc  with  the  moat 
negative  reduced  coat  ia  then  selected  to  enter  the  baaia.  If  no  arc 
has  a  negative  reduced  coat,  then  a  nev  page  (of  area)  is  used.  The 
page  alze  (number  of  area  in  the  page)  is  important.  Too  small  a 
page  will  cause  too  many  pivots j  too  large  a  page  will  cause  too  much 
time  for  the  reduced  cost  calculation. 

After  each  pivot,  a  decision  must  be  made  whether  to  use  the  same 
page  of  data  or  to  obtain  a  nev  page.  One  method  for  making  this 
decision  ia  to  provide  a  re-use  factor,  giving  the  maximum  number  of 
tlmea  a  page  can  be  used  before  a  nev  page  must  be  selected. 

5.1.2  Candidate  List  Methods 

If  arcs  from  a  variety  of  nodes  are  examined,  then,  when  one  of 
them  enters  the  basis,  generally,  it  vill  not  effect  the  reduced  cost 
of  many  of  the  other  arcs.  It  is  therefore  possible  to  examine  only 
a  subset  of  arcs,  called  a  candidate  list.  The  first  step  is  to 
create  a  list  of  arcs  vltlt  a  negative  reduced  cost.  The  arc  vlth  the 
most  negative  reduced  cost  is  then  selected  to  enter  the  basis.  The 
reduced  costs  of  the  arcs  in  the  list  are  then  recalculated,  and  the 
next  arc  to  enter  the  baaia  is  selected  from  the  list.  After  a  fixed 
number  of  pivots,  the  candidate  list  ia  reformed. 

Tvo  parameters  are  required!  the  candidate  list  size  and  the 
number  of  iterations  before  reforming  the  list. 


5. 2  Storing  Arc  Data 


In  both  aajor  nathoda  for  chooaing  an  antaring  arc,  only  a  small 
amount  of  arc  data  is  raquirad  at  any  givan  tiaa.  Information  on 
arcs  in  tha  basis  is  always  raquirad,  but  only  a  small  nuntbar  of  non 
basic  arcs  a ra  naadad.  In  tha  fixad  paga  mathod,  only  arcs  in  tha 
paga  baing  examined  ara  naadad.  In  tha  candidata  list  mathod,  only 
arcs  in  the  candidata  list  ara  raquirad. 

This  suggests  maintaining  arc  data  on  a  high  spaed  mass  storage 
device  (a. g.  hard  disk  on  a  microcomputer > .  Only  basic  arc 
Information  and  a  small  number  of  search  pages  is  maintained  in  core 
When  nev  arc  data  is  required  a  paga  of  arcs  can  be  read  in, 
replacing  the  previous  pages.  Tha  amount  of  data  remains  constant 
while  tha  actual  data  la  constantly  changing.  This  is  called  the 
•in-core/out-of-core"  method. 

Some  method  is  needed  to  store  information  on  tha  non -basic 
flows.  Tha  non-baslc  flows  are  either  zero  or  tha  arc's  upper  bound 
Tha  various  possibilities  ware  given  in  03.  If  there  is  a  large 
number  of  arcs  then  the  information  on  tha  non-basic  flows  must  also 
be  stored  outside  central  memory.  This  is  slow,  but  it  permits  even 
small  computers  to  solve  extremely  large  problems.  * 

9. 3  Specialization  for  HRMATE 

Tha  number  of  arcs  in  a  HRMATE  modal  can  be  vary  large.  Tha 
largest  problems  can  have  more  than  one-half  of  a  million  arcs. 
Problems  of  this  size  require  tha  in-core/out-of -core  method. 
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On*  advantag*  of  th*  HRHATE  probl**  la  that  th*  atructur*  of  th* 
arc  ooata  la  known.  Th*r*  will  b*  aany  area  with  z*ro  coat.  Th*a* 
area  ar*  likely  to  b*  in  th*  optimal  baaia.  Therefor*,  it  aeema 
reaaonabl*  to  enter  zero  coat  area  aa  often  aa  poaalble. 

Th*  remaining  area  can  be  aeparated  into  two  claaaeai  low  coat 
and  high  coat  area.  By  aeparating  th*  area  into  three  different 
file*  (zero,  low,  and  high  coat  filea)  area  with  zero  coat  can  be 
preferentially  entered,  without  calculating  reduced  coata  for  low  and 
high  coat  area.  To  enaure  optimality  all  the  area  muat  be  examined; 
however,  more  time  can  be  apent  with  the  zero  coat  area. 

Slack  and  aurplua  area  ar*  alao  very  Important  in  th*  eolution 
proceaa.  Them*  area  ahould  be  examined  more  often  than  other  area. 

If  th*  "wrong”  alack  and  aurplua  area  ar*  in  the  baaia  many  pivota 
might  be  performed  unneceaaarily.  These  area  ahould  not  be  kept  out 
of  cor*.  Information  necessary  to  generate  these  area  ahould  be 
available  in  core  and  their  reduced  costs  ahould  be  recalculated 
frequently. 


* 
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6.  Finding  th*  Exiting  Arc 

In  th*  simplex  method  speclslized  for  generalized  networks,  only 
a  limited  number  of  arcs  are  candidates  to  exit  the  basis  during  any 
pivot.  The  rapid  identification  of  these  arcs  is  a  reason  the 
specialized  method  is  more  efficient  than  the  simplex  method  for 
general  linear  programming. 

For  any  node,  NODE,  define  the  "backpath*  of  NODE  to  be  those 
arcs  between  NODE  and  the  cycle  for  the  component  that  contains  NODE, 
as  well  as  those  arcs  on  th*  cycle.  In  other  words,  th*  backpath  for 
NODE  contains  those  basic  arcs  whose  corresponding  node  can  be 
reached  from  NODE  by  use  of  the  PRED  structure  only.  For  the 
entering  arc  number  ARC,  the  only  arcs  to  change  flow  are  those  on 
the  backpaths  of  TAIL (ARC)  and  HEAD (ARC).  Figure  6-1  shows  th*  arcs 
that  will  change  flow  for  various  combinations  of  HEAD (ARC)  and 
TAIL (ARC),  referred  to  as  HEAD  and  TAIL  respectively. 


If  a  single  unit  of  flow  were  placed  on  the  arc  entering  the 
basis,  th*  flows  on  th*  arcs  in  the  backpaths  are  the  only  ones  which 
must  be  updated  so  as  to  keep  th*  net  flow  at  each  node  the  same. 

For  instance,  the  arc  between  TAIL  and  PRED (TAIL)  must  provide  on* 
unit  of  flow  at  nod*  TAIL.  If  that  arc  is  orJfcnted  from  PRED (TAIL) 
to  TAIL  and  has  multiplier  of  MULT,  then  the  arc  between 
PRED ( PRED ( TAIL ) )  and  PRED (TAIL)  must  provide  1/HULT  units  of  flow  at 
PRED (TAIL),  and  so  on  (see  Figure  6-1). 

This  calculation  is  equivalent  to  determining  the  updated  column 
in  linear  programming.  With  this  updated  column,  it  is  possible  to 
determine  the  amount  of  flow  by  which  th*  entering  arc  flow  can 
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HEAD 
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Figure  6-1,  Flow  Required  at  Nodes 


ehang*  btlort  mourn  basic  arc  raachaa  on*  of  ita  bounds.  Th*  first 
arc  to  r*ach  on*  of  its  bounds  is  th*  exiting  arc  from  th*  basis. 

6. 1  Calculating  th*  Exiting  Arc 

Th*  d*t*rsination  of  th*  exiting  arc  lnvolvas  two  st*pas  th* 
d*v*lopa*nt  of  th*  updated  column  and  th*  determination  of  th* 
maximum  flow  change  on  th*  entering  arc,  ARC  (before  a  basic  arc 
reaches  on*  of  its  bounds).  Th*s*  two  steps  can  b*  performed 
simultaneously  for  each  nod*.  For  simplification,  th*  presentation 
h*r*  will  separate  th*  process.  There  are  several  cases  to  consider. 

Th*  first  case  assumes  that  the  backpaths  from  TAIL  and  HEAD  are 
distinct.  Implying  that  TAIL  and  HEAD  are  in  different  components. 

If  the  multiplier  on  ARC  is  HULT  then  placing  one  unit  of  flow  on  ARC 
will  require  on*  unit  of  flow  at  TAIL  and  will  provide  HULT  units  of 
flow  at  HEAD.  The  sign  convention  of  a  negative  number  for  a  demand 
and  a  positive  number  for  a  supply  will  be  adopted.  Th*  value  of  >1 
for  TAIL  and  HULT  for  HEAD  is  called  the  requirement. 

Given  th*  requirement  at  a  nod*  HODE,  two  pieces  of  information 
are  required i  th*  update  column  entry  for  the  arc  between  HODE  and 
PRED(HODE)  r.nd  th*  requirement  for  PRED(HODE). 

These  values  depend  on  th#  orientation  of  ARC (HODE),  the  basic 
arc  between  NODE  and  PRED(NODE).  The  following  algorithm  calculates 


th*  entry  in  the  updated  column  (UP_COL)  and  updates  the  requirement 
for  NODE  (REQUIRE)  to  be  th*  requirement  for  PRED(NODE) . 


if  ARC (NODE)  <  0  then  («arc  oriented  from  PRED(NODE)  to  NODE  •) 
REQUIRE  i •  REQUIRE/HULT ( NODE ) ) 

UP_COL ( NODE )  t ■  REQUIRE) 

else 

UP_COL ( NODE )  *»  -REQUIRE) 
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REQUIRE  t-  REQUIRE  •  MULT (NODE); 


endif 

The  proctaf  la  slightly  nor*  complicated  when  a  cycle,  not  a 
sell-loop,  is  reached.  The  effect  of  the  cycle  multiplier  must  be 
taken  into  effect.  Briefly,  the  cycle  is  used  to  create  or  destroy 
flow  as  needed.  The  CMULT  value  determines  the  rate  at  which  flow 
can  be  created  or  destroyed.  The  effect  of  the  cycle  multiplier  is 
that  when  the  cycle  is  reached,  REQUIRE  is  replaced  by  REQUIRE  /  ( 1 - 
CMULT (CYCLE) > .  The  equations  above  can  then  be  used  for  the  arcs  on 
the  cycle.  These  calculations  can  be  carried  out  independently  for 
the  backpaths  of  HEAD  and  TAIL  when  the  nodes  are  in  different 
components. 

For  the  second  case,  when  the  two  nodes  are  in  the  same 
component,  the  backpaths  will  coincide  at  some  point.  The  updated 
column  is  the  sum  of  the  updated  columns  calculated  using  the  above 
equations.  If  the  backpaths  coincide  before  the  root  cycle  it  is 
possible  to  simply  add  together  the  REQUIRE  values  of  the  two 
backpaths  and  continue  as  though  only  one  backpath  existed.  If  the 
value  of  REQUIRE  is  zero  (as  it  will  be  for  pure  networks)  no  further 
calculations  need  be  done;  the  rest  of  the  arcs  will  not  change  in 
flow.  If  the  backpaths  coincide  only  on  the  cycle,  it  is  easiest  to 
proceed  around  the  cycle  twice  making  the  necessary  calculations  and 
add  together  the  resulting  updated  columns. 

An  example  of  these  calculations  is  given  in  Figure  6-2. 

Given  UP_COL(NODE)  it  is  possible  to  determine  the  amount  of 
change  permitted  on  the  entering  arc  before  the  arc  associated  with 
NODE  reaches  a  bound.  Let  INCREASE  i*  1  if  the  entering  arc  is 
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currently  at  zero  flow  end  -1  if  the  entering  arc  hae  flow  equal  to 
lta  capacity.  The  calculation  of  the  maximum  change  la  aa  follow: 


If  (INCREASE  •  UP_COL < NODE )  >  0)  then 

<•  flov  will  decreaae  on  ARC (NODE)  •  ) 

MAX  i *  FLOW (NODE)  /  UP_C0L (NODE > ; 

elee 

(•  flov  vill  lncreaae  on  ARC (NODE)  •  ) 

MAX  *-  (CAPACITY (NODE)  -  FLOW (NODE))  /  UP_C0L(N0DE) » 

endlf 


By  taking  the  alnlaua  value  for  MAX  over  all  nodea  vlth  a 
changing  flov,  the  flov  change  on  the  entering  arc  la  determined  and 
the  exiting  arc  la  identified.  If  thia  value  la  more  than  the 
CAPACITY  of  the  entering  arc  then  the  change  in  flov  ia  the  CAPACITY 
of  the  entering  arc  and  the  exiting  arc  la  the  entering  arc. 
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7.  Updating  the  Baaia 


Glvan  tha  arc  to  antar  and  and  tha  arc  to  laava  tha  baaia,  tha 


final  atap  ia  to  updata  tha  baaia  atructuraa.  Thia  updata  raquiraa 


updating  tha  DUALa,  tha  FLOWa  and  tha  raat  of  tha  baaia  atructuraa. 


7. 1  Updating  tha  Baaia  Structuraa 


7.1.1  Pivot  Typaa 


Tha  linked  rooted  tree  aathod  haa  aix  pivot  typaa,  depending  on 


tha  ralationahip  batvaan  tha  entering  and  exiting  arc. 


Pivot  type  1  occura  vhan  tha  entering  arc  and  exiting  arc  are  tha 


aame.  Thia  occura  vhan  flov  on  tha  entering  arc  raachaa  ita  upper 


bound  or  zero  before  any  other  flova  reach  their  Unite.  In  thia 


caaa,  tha  baaia  ranaina  tha  aame,  ao  only  FLOWa  auat  be  changed. 


Whan  TAIL  and  HEAD  of  tha  entering  arc  are  in  tha  aana  tree,  tha 


pivot  type  la  defined  to  be  either  2,  3  or  4.  Conaidar  patha  from 


tha  entering  tall  node  to  tha  cycle  and  fron  tha  entering  head  node 


to  tha  cycle.  Tha  firat  node  that  occura  on  both  patha  ia  called  tha 


meeting  node  (MEET).  Tha  (common)  cycle  node  ia  called  CYC.  Tha 


exiting  arc  can  occur  in  three  placaai  before  MEET,  batvaan  HEET  and 


CYC,  and  after  CYC.  Thaaa  three  placaa  corraapond  to  pivot  typaa  2, 


3  and  4  raapactivaly. 


Whan  TAIL  and  HEAD  of  tha  entering  arc  are  in  tha  aame  component. 


tha  pivot  type  ia  2  or  S.  If  tha  exiting  arc  ia  on  tha  cycle,  tha 


pivot  type  la  S.  Otharviaa  it  la  type  2. 


If  TAIL  and.  HEAD  of  the  entering  ere  are  in  different  components. 


the  pivot  type  ia  6  if  the  exiting  arc  ia  on  a  cycle,  and  ia  2 
otherwise. 

Figure  7-1  gives  examples  of  all  of  these  pivot  types. 


7.1.2  Common  Routines 

The  sain  advantage  of  the  linked  rooted  tree  method  of  storing 
the  generalised  network  basis  is  that  the  basis  update  routines 
involve  a  limited  number  of  tree  manipulation  routines.  Each  pivot 
type  uses  these  routines  in  a  different  way  to  create  the  new  basis. 

There  are  five  tree  manipulation  routines  required! 

1)  Hang  tree  (HANG)  t  Takes  two  trees  and  a  node  within  the  first 
tree  and  attaches  the  second  tree  to  the  first  below  the  node. 

2)  Isolate  subtree  (ISOLATE)!  Takes  a  tree  and  a  node  within  the 
tree  and  Isolates  the  subtree  below  the  node  by  creating  a  new  tree. 

3)  Reroot  tree  (REROOT)  i  Takes  a  tree  and  a  node  within  the 
tree  and  aakes  that  node  the  root  of  the  tree. 

4)  Reverse  cycle  reroot  (REV_CYC_REROOT) i  Takes  a  series  of  trees 
connected  by  PRED  values  and  creates  a  new  tree  consisting  of  all  of 
them  in  reverse  order. 

5)  Cycle  reroot  (CYC_REROOT) i  Takes  a  series  of  trees  connected 
by  PRED  values  and  creates  a  new  tree  consisting  of  all  of  them. 

7. 1.2. 1  HANG  Routine 


INPUT! 


ROOT:  a  nod*,  not  necessarily  the  root  of  on*  tree; 

TREE:  the  root  node  of  another  tree. 

OUTPUT: 

Updated  PRED,  THREAD,  LEVEL  and  RTHREAD  for  a  tree  vith  the  nodes 


P: 

'  .s 


K  ■- 


f:- 


t  . 

*  L 


toy*# 


of  TREE  belov  ROOT.  <See  Figure  7-2). 


METHOD: 


(A)  Find  the  final  node  In  preorder  traversal  (LAST)  In  the 
subtree  belov  ROOT  by  following  the  THREAD  values) 

(B)  LEVEL (TREE)  :-  LEVEL (ROOT)  ♦  1; 

NODE  : *  THREAD ( TREE ) ) 

while  NODE  <>  TREE  do 

LEVEL (NODE)  :»  LEVEL ( PRED ( NODE ) )  ♦  1) 

NODE  :»  THREAD ( NODE ) ) 
endvhlle; 

(C)  TEMP  :«  THREAD ( LAST ) ; 

THREAD (LAST)  :«  TREE) 

THREAD ( RTHREAD ( TREE ) )  i »  TEMP ) 

RTHREAD (TEMP)  :»  RTHREAD < TREE ) ) 

RTHREAD (TREE)  :«  LAST) 

PRED (TREE)  :»  ROOT) 


DISCUSSION: 

Two  common  programming  structures  are  exhibited  in  this  routine. 
In  section  B,  the  LEVEL  structure  is  updated  for  TREE.  Because  the 
THREAD  is  the  preorder  traversal,  the  level  of  PRED (NODE)  is  always 
calculated  before  the  level  of  NODE.  Since  LEVEL (NODE)  is  always 
equal  to  LEVEL ( PRED ( NODE ) )  ♦  1  within  trees,  the  level  calculation  is 
simplified. 

The  second  structure  is  in  section  C.  For  roots  of  trees,  the 
reverse  thread  of  the  root  is  always  the  last  node  £n  the  preorder 
traversal.  This  means  that  finding  the  last  node  in  TREE  is  easy,  as 
opposed  to  the  last  nod*  belov  ROOT.  Section  A  is  needed  to  find  the 
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Figure  7-2.  HANG  Foutine 
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last  node  below  ROOT 


7.1.2. 2  ISOLATE  Routine 

INPUT i 

NEW_ROOT*  a  node  in  a  tree 
OUTPUT* 

Updated  THREAD  and  RTHREAD  go  that  nodes  below  NEW_ROOT  fora 
a  tree  rooted  at  NEW_ROOT.  (See  Figure  7-3). 

METHOD* 

(A)  Find  the  last  node  (LAST)  in  the  subtree  rooted  at  NEW_ROOT 
by  following  the  THREAD  values} 

(B)  THREAD ( RTHREAD ( NEW_ROOT ) )  *«  THREAD ( LAST ) } 

RTHREAD (LAST)  *-  RTHREAD ( NEW— ROOT ) } 

THREAD (LAST)  *«  NEW  ROOT} 

RTHREAD (NEW_ROOT)  *-  LAST} 

DISCUSSION : 

Once  the  last  node  in  the  subtree  below  NEW_ROOT  has  been 
located,  rethreading  involves  only  nodes  LAST  and  NEW.ROOT. 


BEFORE  ISOLATE 


AFTER  ISOLATE  (1) 


insure  7-3.  Isolate  Routine 


7. 1.2. 3  REROOT  Routine 


INPUT* 

TREE*  the  root  node  of  the  subtree  to  be  rerooted) 
NEW_ROOT*  •  node  in  the  tree  rooted  at  TREE. 

OUTPUT* 

Updated  PRED,  THREAD,  LEVEL  and  RTHREAD  for  a  tree 
consisting  of  the  nodes  belov  TREE,  with  root,  NEW-ROOT.  (See 
Figure  7-4). 


METHOD* 

LEVEL ( NEW  ROOT )  *■  0) 

NODE  *«  NEW  ROOT) 

PREV  *-  NEW  ROOT) 
while  PREV  <>  TREE  do 
ISOLATE (NODE) ) 

TEMP  *-  PRED (NODE)) 
if  NODE  <>  NEW_ROOT  then 
HANG (PREV, NODE) ) 
endif ; 

PREV  *«  NODE) 

NODE  *-  TEMP) 
endwhile) 


DISCUSSION* 

Using  HANG  and  ISOLATE,  rerooting  a  tree  at  a  new  node  is  a 
simple  task.  Each  node  on  the  path  from  NEW_ROOT  to  TREE  is  isolated 
and  then  hung  from  the  previous  node.  The  temporary  variable  TEMP  is 
required  because  HANG  changes  the  value  of  PRED (NODE). 


7. 1.2.4  REVCYCREROOT  Routine 


INPUTS  t 

FIRST_TREE*  Firet  tree  of  the  sequence  to  put  together ; 

LAST_TREE t  Last  tree  of  the  sequence; 

OUTPUTS  t 

Updated  PRED,  THREAD,  RTHREAD  and  LEVEL  values  for  a  single  tree 
rooted  at  FIRST_TREE,  consisting  of  trees  rooted  at  predecessor 
values  fros  FIRST.TREE  to  LAST.TREE.  Note  that  this  definition 
requires  that  the  predecessors  along  the  path  fros  FIRST_TREE  to 
LAST_TR£E  be  reversed.  (See  Figure  7-S). 


METHOD* 


NODE  i-  FIRSTJTREE; 

PREV  *»  FIRST_TREE; 
while  PREV  <>  LAST_TREE  do 
TEMP  *«  PRED ( NODE ) ; 
if  NODE  <>  FIRST.TREE  do 
HANG (PREV, NODE) ; 

endlf 

PREV  *«  NODE; 

NODE  *■  TEMP; 
endvhlle; 


DISCUSSION* 

This  routine  is  essentially  the  same  as  REROOT,  except  that 


ISOLATE  has  already  been  accomplished 


7.  1.2.  S  CYC.REROOT 


:S  i: 


5  8 


a)  Xs 


INPUTS t 


FIRST.TREEx  Firat  t ree  of  the  icqucnc*  to  put  together) 
LAST.TREEi  Laat  tree  of  the  aequence) 


OUTPUTS » 


Updated  PRED,  THREAD,  RTHREAD  and  LEVEL  valuee  for  a  single  tree, 
rooted  at  LAST.TREE,  conaiatlng  of  treea  rooted  at  predeceaaor  valuea 
froa  FIRST.TREE  to  LAST_TREE.  Thla  dlffera  from  REV.CYC.REROOT  only 
In  that  LAST.TREE  la  the  root  lnatead  of  FIRST.TREE.  (See  Figure  7- 


'  S 

f 

: 

* 

3  :•> 


METHOD i 


the  PRED  valuea  from  FIRST_TREE  to  LAST.TREE; 


CYC.REROOT < LAST.TREE, FIRST.TREE) ; 


DISCUSSION  t 


Since  CYC.REROOT  reveraea  the  PRED  valuea,  then  by  reveraing  the 
PRED  valuea  before  calling  CYC.REROOT,  the  PRED  valuea  remain 
unaffected. 


7.1.3  The  Pivot  Routine* 


7. 1.3. 1  Pivot  Type  1 

Since  the  exiting  and  entering  arc  are  the  ■ane,  no  ba*ia 
structure*  other  than  FLOW  must  be  updated. 

7. 1.3.  2  Pivot  Type  2 

Pivot  type  2  can  occur  in  one  of  three  caeeet  entering  head 
and  tail  are  in  the  same  tree,  they  are  in  different  tree*  but  the 
•**e  component,  or  they  are  in  different  component*.  In  all  case*, 
the  exiting  arc  i*  on  either  the  path  from  the  entering  tail  node  to 
the  cycle  or  the  path  from  the  entering  head  node  to  the  cycle,  but 
not  both.  The  LEAVE  node  is  defined  to  be  the  node  associated  vith 
the  path  containing  the  exiting  arc.  The  other  node  is  called  the 
STAY  node  (see  Figure  7-7).  EXIT  refer*  to  the  node  that  i* 
associated  with  the  exiting  arc. 

The  routine  to  execute  a  type  2  pivot  1st 

REROOT (EXIT, LEAVE) } 

HANG (STAY, LEAVE) » 

See  Figure  7-8. 
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Before  Pivot 


KEY 


_  _  _  exiting  arc 
.  .  .  entering  arc 
predecessor 


Figure  7-7.  Pivot  Type  2  -  Initial  Position 
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7. 1.3.  3  Pivot  Type  3 


Pivot  type  3  occurs  vhen  ths  entering  tall  and  head  are  in  the 
ease  tree  and  the  exiting  arc  Is  between  MEET  and  CYC  (Figure  7-9). 
Denote  the  entering  tail  and  head  nodes  as  ETAIL  and  EHEAD. 

The  algorithm  for  pivot  type  3  1st 


NODE  t-  ETAIL; 

PREV  t-  EHEAD; 
while  NODE  <>  MEET  do 
ISOLATE (NODE) ; 

TEMP  »■  PRED(NODE) ; 
PRED ( NODE )  *»  PREV; 
PREV  t-  NODE; 

NODE  t-  TEMP; 
endwhlle; 

NODE  !«  EHEAD; 

PREV  «-  ETAIL; 
while  NODE  <>  MEET  do 
ISOLATE (NODE) ; 

NODE  i-  PRED (NODE); 
endwhlle; 

PRED (MEET)  t-  PREV; 
REROOT ( EXIT, MEET ) 


See  Figure  7-10 
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EHEAD 


Fleure  7-«5 


plvot  Type  3  -  Initial  Position 


/V 


After  Reroot  (5,6) 


Figure  7-10.  *>ivot  Type  3  During  Pivot  (Cont.) 
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7. 1.3.  4  Pivot  Type  4 


Pivot  type  4  occurs  when  the  entering  head  end  tail  nodes  occur 
in  the  same  tree  and  the  exiting  arc  is  on  the  cycle  (Figure  7-11). 
The  algorithm  for  pivot  type  4  1st 


NODE  s«  ETAIL} 

PREY  t»  EHEAD} 
while  NODE  <>  MEET  do 
ISOLATE (NODE) } 

TEMP  »■  PRED(NODE) ; 
PRED ( NODE )  *«  PREY; 
PREY  t-  NODE} 

NODE  t-  TEMP} 
endwhile; 

NODE  t-  EHEAD} 

PREY  *-  ETAIL} 
while  NODE  <>  MEET  do 
ISOLATE (NODE) } 

PREY  t-  NODE} 

NODE  i»  PRED (NODE)} 
endwhile} 

PRED_EXIT  t  *  PRED (EXIT)} 
REROOT (CYC, MEET) } 
REV_CYC_REROOT ( CYC, EXIT) } 
CYC_REROOT(PRED_EXIT, CYC) } 
PRED (MEET)  j-  PREY} 


See  Figure  7-12 
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After  Entire  Routine 


Figure  7-12. 


Pivot  Type  4  (Cont.) 


70 


7. 1.3.  S  Pivot  Type  3 


Pivot  type  3  occur*  when  the  entering  heed  end  teil  ere  in  the 
seme  component,  but  in  different  treee.  The  exiting  ere  occur*  on 
the  cycle  (otherwise  it  would  be  pivot  type  2>.  If  the  predecessor 
psth  from  the  entering  tell  node  to  (end  eround)  the  cycle  is 
compered  to  the  peth  from  the  entering  heed,  one  of  theses  peth  will 
reech  EXIT  first.  Denote  the  node  with  the  peth  to  first  reech  EXIT 
es  LEAVE  end  the  other  es  STAY.  Denote  the  corresponding  cycle  nodes 
ss  CYCJLEAVE  end  CYC.STAY  (see  Figure  7-13). 

The  elgorithm  for  pivot  type  5  is: 


NODE  t-  LEAVE; 

PREV  *-  STAY; 

while  NODE  <>  CYC  LEAVE  do 
ISOLATE (NODE! ; 

TEMP  i •  PRED ( NODE ) ; 

PRED ( NODE >  t«  PREV; 

PREV  i*  NODE; 

NODE  i»  TEMP; 
endwhile; 

PREDJEXIT  *«  PRED (EXIT) ; 
REV_CYC_REROOT ( CYC_LEA VE, EXIT) ; 

NODE  t*  STAY; 

PREV  »«  LEAVE; 
while  NODE  <>  STAY.CYC  do 
ISOLATE (NODE) ; 

PREV  t»  NODE; 

NODE  t*  PRED (NODE); 
endwhile; 

CYC.REROOT ( PRED_EXIT, CYC ) ; 


See  Figure  7-14 


After  Rev-Cyc -Reroot  (3,1) 


Figure  7-14.  wivot  Type  5  During  Pivot 


7. 1.3. 6  Pivot  Type  6 


Pivot  typ*  6  occurs  vh*n  the  entering  arc  has  its  nodes  in 
different  components  and  the  exiting  arc  occurs  on  a  cycle.  <If  the 
exiting  arc  is  not  on  a  cycle  then  the  pivot  type  is  2). 

Denote  the  node  of  the  entering  arc  that  is  in  the  same  component 
as  the  exiting  arc  as  LEAVE  and  the  other  end  as  STAY.  Let  CYC  be 
the  cycle  node  associated  vith  LEAVE.  See  Figure  7- IS. 

The  algorithm  for  pivot  type  6  1st 

REROOT (CYC, LEAVE) ; 

PRED_EXIT  *-  PRED<  EXIT ) ; 

REV  CYC  REROOT (CYC,  EXIT) j 

CYC_REROOT  <  PRED_EXIT, CYC) ; 

HANG (STAY, LEAVE) ; 

See  Figure  7-16. 

7. 2  Updating  the  Duals 

In  order  to  calculate  reduced  costs  quickly,  dual  values 
associated  vith  the  node  are  maintained  in  core.  Because  the  dual 
variable  calculation  is  a  computationally  expensive  operation,  it  is 
fortunate  that  only  a  small  number  of  dual  values  change  at  each 
iteration.  The  duals  that  change  are  exactly  those  whose  node 
receives  a  new  LEVEL  value. 

The  key  to  calculating  duals  is  that  reduced  cost  for 


After  Cyc-Reroot  (4,2) 


After  Pivot 


Figure  7-16. 


Pivot  Type  6  -  During  Pivot 


(Cont , ) 
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arcs  in  tha  basis  are  always  zero.  The  reduced  cost  calculation  is 
<lroa  Section  S) 

MULT  •  DUAL (HEAD)  -  DUAL (TAIL)  -  COST 

Duals  lor  generalized  networks  are  uniquely  determined. 
Calculation  of  duals  for  nodes  on  the  cycle  is  involved  and  is 
reviewed  in  Section  7.2.1.  The  dual  for  a  node  off  the  cycle  is 
baaed  solely  on  the  dual  value  of  its  predecessor.  This  is  shown  in 
Section  7.  2.  2. 


7. 2. 1  Dual  Values  on  the  Cycle 


Given  a  cycle  with  k  nodes  there  are  k  arcs  between  them.  Each 
arc  creates  one  reduced  cost  calculation.  Therefore,  there  are  k 
linear  equations  to  find  k  unknown  duals. 

If  one  dual  on  the  cycle  is  given,  the  rest  can  be  obtained 
by  traversing  the  cycle  with  the  PRED  values.  The  method  to  find 
one  dual  on  the  cycle  follows. 

Let  the  nodes  on  the  cycle  to  be  1  . . .  k,  with  PRED(i)  *  i ♦ 1  and 
PRED(k)  *  1.  The  following  calculates  DUAL(l): 


TOT  SUM  i«  0; 

TOT_HULT  i«  If 
NODE  f  1 ; 

repeat 

if  ARC (NODE)  <  0  then  (*  ARC  is  from  PRED (NODE)  to  NODE  •) 
TOT_MULT  »■  TOT  MULT  /  MULT ( NODE )| 

TOT_SUM  f  TOT_SUM  -  TOT_MULT  •  COST (NODE); 

else 

TOT_SUM  *«  TOT_SUM  ♦  TOT_MULT  *  COST (NODE); 

TOT_MULT  t-  TOT.MULT  •  MULT (NODE); 
endif ; 

NODE  t-  PRED (NODE); 
until  (NODE  «  1); 

DUAL ( 1 )  «■  TOT_SUM  /  ( 1-CMULT( 1 ) ) ; 


l*r 


s 


I 
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Duals  of  cycle  nodes  associated  with  self-loops  are  easier.  For 
this  case,  the  dual  Is  simply  the  cost  divided  by  the  self-loop 
multiplier. 

7. 2. 2  Dual  Values  not  on  the  Cycle 

The  dual  value  for  a  node  off  the  cycle  is  calculated  asi 

if  ARC (NODE)  <  0  then 

DUAL (NODE)  f  ( DUAL (PRED< NODE ) )  ♦  COST)  /  MULT; 

else 

DUAL (NODE)  i»  DUAL (PRED( NODE ) )  *  MULT  -  COST; 

endlf 

It  is  important  that  the  dual  variable  updates  are  performed  in 
the  correct  order.  For  every  node,  the  dual  for  the  PRED  of  the  node 
must  be  calculated  before  the  dual  of  the  node  can  be  calculated.  It 
is  for  this  reason  that  THREAD  is  defined  to  be  the  preorder 
traversal.  Following  the  thread  will  ensure  update  of  the  PRED  of 
any  node  before  the  node  itself ;  the  DUAL  update  can  be  accomplished 
at  the  same  time.  This  reduces  the  number  of  times  any  node  must  be 
examined  during  a  pivot,  yielding  in  a  significant  reduction  in 
computation  time. 

7. 2  Updating  the  Flows 


The  final  structure  to  update  is  the  basic  flow  values.  Host  of 
the  work  in  this  update  was  accomplished  during  the  calculation  of 
the  exiting  arc.  The  updated  column  was  the  change  in  flow  on  the 
basic  arcs  if  one  unit  of  flow  was  put  on  the  entering  arc.  This 
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allowed  a  calculation  of  the  maximum  flow  that  could  be  placed  on  the 
entering  arc  without  violating  any  bounda.  Thla  value  la  multiplied 
by  the  updated  column  and  added  to  the  current  flowa  to  produce  the 
updated  flowa. 

Algorithmically,  It  la  eaaier  to  perform  thla  update  while  the 
baaia  la  being  updated.  A a  each  node  la  vialted  It  la  a  simple 
natter  to  calculate  the  amount  the  flow  will  change  and  update  FLOW 
accordingly.  Some  nodea  are  not  vialted  during  the  baaia  update  and 
thoae  nodea  nuat  be  vialted  aolely  to  update  the  flowa.  For 
inatance,  in  pivot  type  1,  no  baaia  update  la  performed  but  the  flowa 
muat  be  updated.  All  nodea  vialted  in  the  calculation  of  the  exiting 
arc  muat  be  revlalted  to  update  the  flowa. 


8.0  OTHER  CONCERNS 


Major  portions  of  tha  primal  aimplax  aathod  hava  baan  covarad. 
Givan  a  basis  it  is  possibla  to  find  an  arc  to  antar  tha  basis, 
datarmina  tha  arc  to  axit  tha  basis  and  updata  tha  basis.  Tha  only 
othar  stap  of  tha  primal  aimplax  mathod  is  to  find  an  initial  basis. 
Saction  8. 1  gives  tvo  aathods  for  doing  this. 

Ona  othar  concarn  is  tha  afficiancy  of  tha  generalized  natvork 
primal  aimplax  aathod  vhen  solving  a  pura  natvork.  It  is  genarally 
thought  that  a  pura  natvork  coda  vill  axacuta  tvo  or  thraa  timas 
faster  than  a  ganaralizad  natvork  coda.  Soma  of  this  improvement  is 
dua  to  tha  fact  that  ganaralizad  natvorks  must  kaap  track  of  tha 
multipliars  vhile  tha  pura  natvork  has  ona  lass  piaca  of  data  for 
aach  arc.  Some  afficiancy  ralatas  to  tha  structura  of  tha  pivota 
usad  vith  pura  natvorks.  This  pivot  and  basis  structure  can  be  used 
vithin  a  ganaralizad  natvork  coda,  so  pura  natvorks  and  "almost"  pure 
natvorks  can  ba  solvad  mora  quickly  than  generalized  natvorks  vith 
tha  generalized  natvork  simplex  mathod.  This  is  detailed  in  Saction 
8.2. 


8. 1  Initial  Basis 


Tha  starting  basis  has  a  large  effect  on  time  required  by  tha 
primal  aimplax  mathod.  Ideally,  if  tha  optimal  basis  is  usad  as  tha 
starting  basis,  than  no  pivota  need  ba  performed.  Tha  concept  of 
quality  of  tha  starting  basis  is  difficult  to  quantify,  although  it 
seams  reasonable  to  assume  that  a  good  starting  basis  is  ona  vhlch 


has  a  large  number  of  area  that  are  also  In  the  optimal  baala. 

If  too  much  time  la  spent  In  trying  to  find  a  good  starting 
basis,  the  total  execution  time  might  be  larger  than  using  a  poorer, 
but  easier  to  find,  basis. 

An  initial  basis  must  satisfy  the  following  conditions* 

1)  There  must  be  one  arc  in  the  basis  for  each  node. 

2)  Each  component  generated  by  the  arcs  in  the  basis  must  have 
exactly  one  cycle,  which  nay  be  a  self-loop. 

3)  The  net  flows  into  and  out  of  each  node  must  be  equal  to  the 
supply  or  demand  at  that  node. 

4)  The  flow  on  any  arc  not  in  the  basis  must  be  either  zero  or 
the  capacity  of  the  arc. 

5)  The  flow  on  any  arc  must  be  between  zero  and  the  capacity  of 
the  arc  inclusive. 

Section  8.1.1  describes  the  easiest  basis  to  find*  the 
artificial  start  basis.  By  spending  essentially  no  time  in  creating 
the  initial  baala,  this  basis  makes  no  attempt  to  guess  which  arcs 
will  be  in  the  optimal  basis. 

Section  8.1.2  discusses  basis  creation  methods  that  try  to  guess 
which  arcs  are  optimal.  These  advanced  start  methods  are  often 
problem  specific,  although  some  general  purpose  methods  are  possible. 

8. 1. 1  Artificial  Start 

Associated  with  each  node  in  a  generalized  network  is  a  self¬ 
loop,  representing  the  slack  or  surplus  variable  at  the  node.  If  the 
supply  or  demand  at  the  node  must  be  satisfied  then  there  is  a  large 
coat  on  the  self-loop,  otherwise  the  cost  is  zero.  The  artificial 


•tart  baala  la  as  follows* 

1)  The  flov  on  avary  self-loop  la  tha  supply  or  daaand  at  the 
noda. 

2)  Tha  flov  on  avary  othar  arc  is  zaro. 

Tha  area  in  < 1 )  abova  fora  the  baais.  This  basia  aatlaf laa  tha 
f lva  eonditlona  givan  In  Section  8. 1.  It  la  alao  eaay  to  create. 
Furtheraore  the  PRED,  THREAD,  LEVEL  and  RTHREAD  are  very  eaay  to 
calculate.  Figure  8-1  llluatratea  thla  for  tha  network  in  Figure 
2-3. 

Tha  diaadvantage  of  thla  aathod  la  that  no  attaapt  la  aade  to 
dateraina  which  area  will  likely  occur  in  tha  optiaal  baala.  Tha 
nuabar  of  plvota  required  after  an  artificial  atart  la  probably  aore 
than  tha  nuabar  required  after  other  types  of  starts. 

8.1.2  Advanced  Start 

A  baala  that  atteapts  to  contain  arcs  that  will  occur  in  the 
optiaal  baais  la  referred  to  as  an  advanced  start.  Advanced  starts 
can  often  be  deterained  froa  the  structure  of  the  particular  problem 
types  that  are  being  solved.  Many,  although  not  all,  advanced  start 
strategies  adopt  tne  following  form 

1 )  Sort  the  area  in  decreasing  order  of  attractiveness 
(likelihood  to  be  in  optimal  baais) 

2)  Using  the  moat  attractive  unexamined  arc,  place  the  maximum 
amount  of  flov  that  will  keep  feasibility  at  the  nodes  and  the  arc  on 
the  arc. 

3)  Repeat  Step  2  until  all  of  the  arcs  have  been 

examined. 

Algorithms  that  adopt  the  above  strategy  and  place  all 
unallocated  supply  and  demand  on  the  self-loops  create  a  basis  that 


Figure  8-1,  Artificial  Basis 
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satisfy  ths  fivs  conditions  in  Ssction  8.1  (possibly  sons  sslf -loops 
with  zsro  flov  will  havt  to  bs  sddsd  to  ths  basis).  If  sasy  nsthods 
of  calculating  attractivsnssa  ars  availabls  and  if  attractivsnsss  is 
a  trus  nsasurs  of  ths  likslihood  of  bsing  in  ths  optimal  basis  than 
advancsd  starts  of  this  forn  ars  gsnsrally  fairly  sffsctivs. 

In  Figurs  8-2  ths  arcs  in  Figurs  2-3  ars  ordsrsd  by  incrsasing 
cost  and  an  advancsd  start  basis  is  formsd. 

8.1.3  Specialization  for  MRMATE 

Arcs  of  ths  MRMATE  aodel  ars  already  sorted  into  three  classes 
based  on  attractiveness.  This  pernits  an  sasy  advancsd  start.  Arcs 
in  ths  zero  coat  arc  file  are  placed  into  the  basis  first.  Next,  low 
cost  arcs  are  used.  Because  there  are  a  large  number  of  high  cost 
arcs,  it  seems  better  to  ignore  them  in  creating  an  advanced  start 
and  let  the  rest  of  the  simplex  algorithm  choose  which  high  cost  arcs 
to  use  by  pivotting  them  into  the  basis  as  needed. 

8. 2  Effect  of  Pure  Network  Structure 

Pure  networks  can  bs  solved  more  quickly  than  generalized 
networks.  This  is  for  two  reasons < 

1)  Pure  networks  have  multipliers  of  one.  Hence,  rather  than 
using  multiplications  and  divisions,  the  various  simplex  calculations 
can  assume  the  MULT  is  one. 

2)  The  basis  structure  of  pure  networks  is  simpler  than  that  of 
generalized  networks.  Each  component  of  a  pure  network  basis  must 
be  rooted  at  a  self-loop. 
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It  la  possible  lor  a  generalized  network  code  to  solve  pure 
networks  alaost  as  quickly  a  pure  network  code  by  perforalng  the 


following t 

1>A  aultlpller  Is  compared  to  1  before  a  division  or 
multiplication.  If  It  equals  1  then  the  division  or  multiplication 
Is  not  needed. 

2)  Data  structures  are  used  that  are  equivalent  to  pure  network 
data  structures  in  the  case  of  a  self-loop  root. 

Unfortunately,  handling  the  multipliers  as  suggested  In  point  one 
will  slow  down  execution  for  all  networks.  This  decrease  In  speed 
may  not  be  excessive  but  may  not  be  worthwhile  if  the  number  of  pure 
arcs  is  small.  Speed  increases  will  occur,  however,  even  in 
generalized  networks,  as  long  aa  the  network  has  a  large  proportion 
of  pure  arcs. 

Data  structures  used  in  pure  networks  are  exactly  those  outlined 
in  Section  4  in  the  case  where  the  root  cycle  is  a  self- loop. 
Advantages  of  the  basis  structure  can  occur  in  any  generalized 
network  that  can  be  converted  to  a  pure  network  by  scaling  the  rows 
and  columns  of  the  constraint  matrix.  The  basis  advantages  can  be 
gained  even  if  the  scaling  is  not  performed.  An  example  of  a 
generalized  network  that  can  be  scaled  to  a  pure  network  is  given  in 
Figure  8-3. 

The  main  advantage  of  the  pure  network  structure  is  that  all 
pivots  will  be  of  types  1  or  2  (if  the  self -loops  are  replaced  with 
artificial  arcs  and  an  artificial  node,  see  C3J).  No  cycles  are 
ever  created,  so  no  cycle  multipliers  need  be  calculated  (Section  4). 
Calculation  of  the  cycle  multiplier  is  a  very  expensive  operation, 
since  it  Involves  many  multiplications  and  divisions. 
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It  seems  likely  that  a  network  that  ie  "almost"  transformable  to 
a  pure  network  would  have  many  of  the  same  advantages  as  a  network 
transformable  to  a  pure  network.  "Almost"  must  be  defined  very 
carefully.  For  instance.  Figure  8-4  shows  a  network  with  just  one 
arc  with  a  multiplier  not  equal  to  one.  However,  this  network  is 
equivalent  to  the  network  in  Figure  2-3  and  the  two  networks  are 
solved  in  almost  the  same  way  by  a  generalized  network  simplex 
method.  The  number  of  cycle  calculations  is  the  same,  as  is  the  time 
required.  The  reason  for  this  is  that  non-pure  arcs  occur  in  the 
basis  disproportionately  for  their  quantity.  This  shows  the 
disadvantage  of  defining  "almost*  pure  networks  in  terms  of  the 
percentage  pure  arcs. 

"Almost"  pure  networks  are  those  whose  valid  basis  tend  to  have 
self-loops,  rather  than  cycles.  Cycles  are  formed  when  the  exiting 
arc  occurs  on  the  common  part  of  the  backpaths  formed  by  the  entering 
arc.  "Almost"  pure  networks  require  that  for  a  "typical”  basis  and 
entering  arc,  the  exiting  arc  occurs  on  the  separate  parts  of  the 
backpath  (see  Section  6.1). 

8. 2. 1  Specialization  for  HRHATE 

As  might  be  expected,  the  HRHATE  model  can  be  a  pure, 
transformable  to  pure,  or  "almost”  transformable  to  pure  network.  If 
there  are  no  sea  channels  then  all  of  the  multipliers  in  the  network 
are  one,  so  a  pure  network  results.  If  there  are  no  air  channels 
then  the  network  is  transformable  to  a  pure  network.  This  is 
equivalent  to  using  the  volume  of  the  movement  requirement  rather 
than  the  weight. 


If  there  are  far  more  air  channels  than  sea  channels,  or  sore  sea 
channels  than  air  channels,  then  the  resulting  network  is  "almost” 
transformable  to  a  pure  network  as  long  as  the  less  numerous  channels 
are  not  disproportionately  important.  This  is  due  to  the  following. 
Given  a  network  with  far  more  air  channels  than  sea  channels,  a  basis 
and  an  entering  arc,  if  the  channels  are  equally  important,  the 
typical  entering  arc  will  be  an  air  channel  arc.  The  backpaths 
represent  the  channels  that  will  have  to  change  flows  in  order  to 
permit  the  entering  arc  to  Join  the  basis.  Typically,  thlt,  will  only 
involve  air  channels.  This  means  that  the  exiting  arc  cannot  occur 
on  the  common  part  of  the  backpath,  since  a  cycle  with  a  cycle 
multiplier  of  1  would  be  created.  A  similar  argument  holds  for  more 
sea  channels  than  air. 

How  many  more  is  "far  more"?  Certainly  75  air  channels  and  25 
sea  channels  will  not  exhibit  much  pure  network  structure,  since  moat 
backpaths  will  contain  both  air  and  sea  channels.  Also,  99  sea 
channels  and  1  air  channel  will  exhibit  a  lot  of  pure  network 
structure,  since  most  pivots  will  reassign  flow  among  the  sea 
channels.  The  effect  of  pure  network  structure  is  further  examined 
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